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1. Introduction 

The Landau-Lifshitz (LL) model describes the dynamics of the classical Heisenberg spin chain, and has 
appeared in the past few years in both gauge and string theories as one of the first indications of the underlying 
integrable structure [3-15] of the AdS/CFT correspondence (for a review see [16-19]). On the string theory 
side the LL model emerges in the R x subsector in the limit of large angular momentum [3, 20-25]. 
Alternatively, the LL model can be derived from the Faddeev-Reshetikhin (FR) model in the low-energy 
limit [26,27]. The latter is aimed to resolve the difficulties of the quantization procedure [28-36] inherent to 
all sigma models. More recently, the LL model also emerged in the SU{2) x SU{2) subsector of the strings on 
AdSi X CP^ [37] in the context of the newly proposed duality between M = 6 superconformal Chern-Simons 
and the ^4^54 x CP^ string theories in the t'Hooft limit [38]. The integrable properties of the LL model 
have been extensively discussed from various points of view. The classical integrability was established for 
the isotropic case in [39], and for the general anisotropic case in [29,40], and the quasi-classical analysis was 
performed in [41]. The classical equivalence between the LL and the non-linear Schrodinger (NLS) models was 
shown in [42] by constructing a gauge transformation between the ffat currents of the corresponding models. 

At the quantum level, on the other hand, the integrable properties of the LL and NLS models are quite 
different. In [43] Sklyanin has considered in detail the quantum inverse scattering method for the LL model 
(the su{l, 1) case), and pointed out several subtleties of the quantization procedure, absent in the classically 
equivalent NLS model. One of the surprising difficulties is that the standard methods to obtain the Yang- 
Baxter relation from its classical counterpart fail for the anisotropic LL model, and the algebra of observables, 
as well as the transfer matrix have to be modified by hand in order to construct the quantum R-matrix. More 
important implication of the Sklyanin's analysis is the problem of constructing the local integrals of motion in 
the quantum theory, underlying the quantum integrability of the model. In fact, the problems already arise 
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in constructing the local quantum Hamiltonian for the two-particle sector, and to the best of our knowledge, 
the higher order charges have not been constructed so far. The difficulties arise as a result of the ill-defined 
operator product at the same point^ and can, in principle, be resolved by the standard procedure of putting 
the model on the lattice [45,46]. Although this can be done for both FR and LL models [26,47,48], and the 
corresponding spectra can be found with the standard Bethe Ansatz technique, the full AdS^ x string model 
is too complicated, and in general putting a continuous theory on the lattice using the standard methods leads 
to a non-local quantum Hamiltonian [45,49]. Thus, one has to deal with the quantization of the continuous 
models directly, and the LL model, despite its simplicity, is a particularly interesting representative example 
of the difficulties associated with the continuous quantum inverse scattering method. 

The LL and FR models, as well as a number of other models, arising in various limits of the full AdS^ x 
string, have recently been considered as two-dimensional field theories, and the corresponding S-matrices have 
been obtained by means of perturbative calculations [27,50-53]. The idea of this method and its relation to 
the inverse scattering method is nicely reviewed for the NLS model in [54]. This method makes use of the 
integrability, which is believed to hold at the quantum level, and, in particular, the S-matrix factorization 
property [55-62], which renders the N-particle scattering S-matrix to be a product of the two-particle S- 
matrices. The direct verification of the factorization property, by calculating the set of Feynman diagrams, 
has only been possible to carry out for the NLS model in [1,2]. More recently, the three-particle S-matrix 
factorization property was demonstrated for the AdS^ x string in the near fiat space limit at one-loop 
order [53]. 

In this article we consider the three-particle scattering process for the LL model and, by computing the 
necessary Feynman diagrams up to the second order, show that it is factorizable in the first non-trivial order. 
There are several conceptual features that make the S-matrix calculations for the LL model a non-trivial task. 
First, as we noted earlier, the LL model is the low energy limit for the FR model, and as was shown in [36], 
the Hamiltonian for the FR model is not diagonalizable in the class of the standard representation for the 
two-particle sector: 

This representation is at the heart of the Bethe Ansatz and has the clear interpretation for the first and 
second terms as the incoming and outgoing waves respectively, and where the S-matrix S{pi,p2) = e*^*^^^'^^) is 
a phase-shift due to the scattering (we assume that pi > P2)- Proceeding in a similar manner to [36] it is not 
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difficult to show that the Hamiltonian for the LL model cannot be diagonalized for the class of functions (1.1 ), 
although the exact two-particle scattering S-matrix for the FR model can be obtained using the perturbative 
calculations [27]. Alternatively, this can also be seen directly via the quantum inverse scattering method [43], 
where the two-particle wave function is shown to acquire an additional term to (^]^). One of the reasons why 
this happens is the extremely singular nature of the LL interaction [43]. In fact, it is not difficult to see that 
the interaction in the quantum mechanical picture corresponds to the second derivative of the delta-function, 
and dealing with it is not an easy task (see for example [63]). In order to write the additional term to 



(1.1), one has to introduce a new creation operator, which for the two-particle state would correspond to a 
bound state. It is not difficult to see that the number of such terms increases with the number of particles 
in the scattering process. For example, in the three-particle scattering process there will be two additional 
contributions, corresponding to three and two-particle clusters. 

^For more details see the upcoming paper [44] 
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More importantly, the form (|1 . 1| ) suggests that the particles created by the fields 0^(x)|O) are those 
corresponding to the Bethe particles, and all the consequences of the integrability, such as the S-matrix 
factorization, will hold as the consequence of the general argument of [55,56,59]. This is indeed the case for a 
number of models, e.g. the NLS model, where ( |1.1| ) is indeed the two-particle state, and the particles 0"^(x)|O) 
are the Bethe particles, for which one can obtain the S-matrix either using the perturbative diagrammatic 
calculations, or via direct diagonalization of the Hamiltonian. However, this a priori is not the case in general, 
and the particles created by the fields (p~^{x) may not correspond to the Bethe particles. This is the case for 
the LL model, namely, the Bethe particles do not coincide with (x)|0). Let us emphasize that considering 
the scattering process and stating the S-matrix factorization makes sense only in terms of the Bethe particles. 
The only reliable method for diagonalizing all the local conserved quantities simultaneously, including the 
Hamiltonian, and, as a consequence, constructing the eigenstates corresponding to the Bethe particles, is the 
quantum inverse scattering method. But as we mentioned above, this construction fails in the LL model, and 
the local conserved quantities cannot be derived from the trace identities due to the singular expressions. Thus, 
checking the quantum integrability, e.g., the factorization of the S-matrix, as well as establishing a connection 
with the construction in the quantum inverse scattering method by using the standard field theoretic methods 
is an important task that should prove useful when considering more complicated systems than the LL model. 
In this paper and [44] we make the first steps towards this program, and as the first result we show here the 
factorization in the lowest order. 

There are other subtleties, making factorization for the LL model a non-trivial feature. The unconstrained 
two dimensional field theory on the world-sheet corresponding to the LL model contains infinite number of 
vertices. This is in contrast with the simpler NLS model, where there is simply one vertex to deal with, and, 
thus, making the higher order calculations much simpler. In the LL model one has to correspondingly consider 
new types of vertices as the order of calculations increases. Let us recall, that to calculate the two-particle 
scattering S-matrix one only needs to keep the terms up to the quartic order in the Lagrangian. The role 
of higher order terms is to preserve the integrability at the quantum level. Namely, for the three-particle 
scattering S-matrix one has to consider the vertices of higher (up to cj)^) order. The quantum integrability 
manifests here in the S-matrix factorization for the three-particle scattering. Although the combinatorial 
analysis due to the infinite number of vertices makes practically impossible to carry out the calculations to all 
orders, we are able to show in the first non-trivial order that the higher order vertices are exactly of the form 
needed for the S-matrix to be factorizable. The mechanism behind factorization is not straightforward, and we 
show that the higher order terms cancel unwanted lower order terms to guarantee the factorization. Moreover, 
we show explicitly that there is no process corresponding to particle annihilation or creation, and the set of 
momenta before and after the scattering is the same, in agreement with our understanding of integrability. 

Our paper is organized as follows. In Section |^, we give a brief account of the perturbative S-matrix 
calculations for the Landau-Lifshitz model in the two-particle sector. In Section [s], we present our analysis 
of the three-particle scattering for the LL model and show the factorization in the first non-trivial order. In 
section ^, we give a brief summary of our results. In the Appendix we collect the set of the Feynman graphs 
used in our calculations in order to avoid cluttering the main text. 

2. Two-particle scattering S-matrix 

In this section we set up the notation and present the necessary formulas for the subsequent section on the 
three-particles scattering. First we briefly review the necessary facts about the LL model, and then explain 
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the two-particle perturbative S-matrix calculations, referring the reader to [27] for complete details of the 
calculations. Let also us note, that the two-particle scattering diagrammatic calculations for the LL model 
are quite similar to those of the non-linear Schrodinger model [54]. Indeed, the only essential difference is 
the presence of the derivatives in the interaction terms. In the context of the perturbative calculations this 
results in a slightly more complex combinatorial analysis in the LL model. As we will see below, the essential 
differences arise when considering the three-particle scattering amplitude. 

2.1 Landau-Lifshitz model: preliminary facts 

The Landau-Lifshitz model is a non-relativistic sigma model which describes the continuous Heisenberg fer- 
romagnet, and arises in the R x subsector of the AdS^ x strings in the limit of the large angular 
momentum. The equations of motion for the isotropic LL model, on which we will focus in this paper, in 
terms of the spin variables S{x) = (5"(x); a = 1, 2, 3) , have the form: 

dtS = SxdlS (2.1) 

where the fields S(x) take values on S^: 

= 1 (2.2) 
The equation ( |2.1| ) can be obtained from the Hamiltonian: 

HLL = \j dx{d^Sf (2.3) 

with the Poisson structure: 

{S''{x),S\y)] = -e'^'''S'^5{x-y) (2.4) 

In the most general formulation, the equations of motion for the anisotropic LL model are modified by an 
anisotropy matrix J = diag(Ji, J2, Js), and take the form: 

5tS = Sx52S + S X JS (2.5) 

As explained in [43] the isotropic and anisotropic cases should be considered separately when quantizing the 
theory. In the former case, the standard inverse scattering procedure goes through without any changes, and 
the Yang-Baxter equation is satisfied with the appropriate choice of the i?-matrix, while in the later case 



the operators, and accordingly the operator algebra (2.4) have to be modified by hand for the Yang-Baxter 



relation to have a solution. In this paper we will consider only the isotropic case J = 0, and the more general 
case is currently under investigation. 

The corresponding action can be written in a manifestly covariant form as follows: 



S= I d^x 

where Ct(S) is the Wess-Zumino term: 



(2.6) 



Ct{S) = -I C di eijkSid^SjdtSk (2.7) 
^ Jo 

The boundary conditions for the S(t, x;^) field have the form: 
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1 S{t,x;C = 1) = So . . 

\ S{t,x;C = 0) = S{t,x) ^ ■ ^ 

where So is a constant vector. 

Following [27], it is convenient to resolve the constraint ( |2.2D and write the action ( |2.6| ) as an unconstrained 
(1 + l)-dimensional quantum field theory. This can be achieved by a change of variables, which get rid of the 
non-linearities in the kinetic term, as follows 



and which turns the Wess-Zumino term ( |2.8D into a but non-covariant form. The resulting action following 
from ( |2.6D and (|2.9| ) can be written as follows: 



S = ld'x {^(^*9*c,-9,c,V)-|5.v,|2-i^-^[(^*9.^)2 + a^^^ (2.10) 

The non-relativistic character of the LL model leads to significant simplifications in the quantization 
procedure. In particular, it makes diagrammatic calculations quite easy to deal with. This is mostly due to 
the property of the propagator in the LL model to be a retarded one. Let us also note, that in relativistic 
theories, such as the massive Thirring model, to have similar simplifications in diagrammatic calculations, 
one can choose the false vacuum by hand, and, therefore, making the Feynman propagators to be retarded. 

Using the field decomposition 

^{x) = I ^e-P-ap , ^*{x) = I ^e^^-^al (2.11) 

where we use the notation p ■ x = p^x^ — px, and defining the vacuum of the theory by ip{x)\0) = 0, the 
annihilation and creation operators Op and a], satisfy the standard commutation relation: 

[ap,,al]=27T6{p-p'). (2.12) 
One then finds the propagator to be of the form [27]: 



d q ie ^'^'^ 



(27r)^ q^-q^+ie 

(2.13) 



D{t,x) 



The fact that the propagator is purely retarded leads to the following result: the two particle scattering 
S-matrix is given by the sum of bubble diagrams, as in Fig. |l]. 

The diagram above is exactly the one that appears in the NLS model. For the LL model one has to 
properly take into account the presence of the derivatives in the interaction terms, which results in four 
different types of vertices, depending on the placement of the derivatives on the internal or external lines. 
Nevertheless, the calculations are still easy to carry out. 
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P k 
Figure 1: Typical bubble diagram for the two-particle S- matrix. 

2.2 Diagrammatic calculations 

Since the only contributions to the two-particle scattering amplitude are the diagrams of the type Fig. |l|, to 
calculate the two-particle S-matrix it is enough to truncate the complicated action ( p.lO| ) up to the fourth 
order in the fields: 

S = j Sx '-{ip^dtif - dtip\) - \d,ip\^ - |[(^* 5,95)2 + {d,ip*ipf] + 0(996) (2.14) 

We introduced in the expansion ( |2.14| ) the formal parameter g to keep track of the perturbative order, and 
it should be set to one at the end of the calculations. Let us note here, that this is very important for the 
three-particle scattering amplitude factorization, where several such formal parameters will be introduced to 
take care of order counting for different types of vertices that arise in the LL model when expanding the 
action ( |2.10| ) up to the sixth order in the fields. A natural question arises whether setting the parameters to 
one at the end is consistent with the renormalization properties of the LL model. Although we ignore here 
all the divergences associated with the loop calculations, this is an interesting question and is currently under 
investigation. 

We collect here the necessary expressions, referring the reader to [27] for complete details. Taking without 
any loss of generality the ordering p > p' for the scattering to take place, the two-particle S-matrix is 
determined from the relation 

{kk'\S\pp') = S{p,p')8f{p,p';k,k') (2.15) 

where: 

(^i^^ (p, p'- k, k') = (27r)2 [5{p - k)6{p' - k') + 6{p - k')5{p' - k)] (2.16) 
The scattering amplitude given by: 

{kk'\S\pp') = {kk'\Te-''f"^^''^^\pp') 

= {kk'\pp') - i{kk'\T J Hintdt\pp') - ]^{kk'\T [^j Hintdt^ \pp') + • • • (2.17) 

is easily computed with Hint = f / \{'^* ^^'■P)'^ + {^x'^*^)'^'] in each order. The non-scattering part and the 
tree level parts are easily computed and have the forms: 

{kk'\pp') = 6f{p,p';k,k') (2.18) 

-i{kk'\T [ Hintdt\pp') = 2ig^^5^^\p,p';k,k') (2.19) 
J p-p' ^ 
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To avoid the combinatorial analysis of [27], associated with placement of derivatives on internal and external 
lines, we will compute the full bubble diagram in the following manner. Using the complete interaction vertex, 
which in the momentum representation has the form: 

V{k, k';p,p') = 2ig [pp' + kk'] 6{p^ + - k^ - k'^)6{p + p' -k-k'), (2.20) 

we compute the complete one-loop diagram, while keeping the external line corresponding to k and k' off-shell, 
and putting the momenta p and p' on-shell. The resulting expression has the following form: 

{kk'\S\pp')\g., = \2ig {pp' + kk') {27r)'^6^^\p + p' - k - k')] ig^^{pp' + kk') (2.21) 



P-P 



This already includes the sum of all one-loop diagrams with all possible placement of derivatives on internal and 
external lines, with the correct combinatorial factors. Noting that the term inside the square brackets in ( p. 21 ) 



coincides with the expression for the interaction vertex ( p.20| ) , one can regard the one-loop scattering amplitude 
as the interaction vertex in momentum representation, with the momenta p and p' on-shell, multiplied by 
some function of this pair of momenta and the expansion parameter. So that the n-loop scattering amplitude 
corresponds to the product of n of these modified vertices. 

i9j^] 2ig{27Tfd^^\p + p'-k-k'){pp' + kk') (2.22) 
If we also put k and k' on-shell, we find that: 

(/ \ n+l 
^9^) S^^\p,p';k,k'), (2.23) 

where we have used the relation 

(27r)25(2) + y _ ^ _ ^/) ^ 1 .(2) p>. ^/) (2.24) 

2[p — p'j 

Thus, two-particle scattering S-matrix has the form [27]. 

Sip^p') = 1 + 2 f:i^gr = i^ir^ ^^fe (^-^^^ 

1 \P-pJ --^+iq --^+1 

n=l \i- 1- / p p' ' ^ p p' ' 

3. Three-particle scattering S-matrix 

The direct verification of the S-matrix factorization using the perturbative calculations is quite difficult to 
carry out, due to the complicated diagrammatic analysis. The only model for which the factorization property 
was manifestly verified is the non-linear Schrodinger model. We refer to the details of the diagrammatic 
calculations for the case of the NLS model to [1,2]. Before presenting the concrete calculations for the LL 
model let us note that there is a significant difference between the NLS and LL models, and the latter case is 
considerably more complex. Indeed, the Hamiltonian for the NLS model has the form: 

H= [ dx \d^(Pix)\'^ + c\(Pix)\^ (3.1) 
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Thus, while calculating the three-particle scattering amplitude, one has to deal with only one type of vertex. 
The same is true when considering the general A^-particle scattering process. In the LL model the situation is 
different and to consider the three-particle scattering amplitude one has to truncate the full action ( |2.10| ) up 
to the sixth order in fields. The higher order vertices will not contribute as the consequence of the retarded 
propagator as well as the charge conservation [27]. The resulting Lagrangian has the form 



C='-i^*dtV-dtip*v)-\dM'-Y^i^'^.^)' + (^-^*^)'^-^\^^^^^^ (3-2) 

where we have introduced, analogous to the two-particle scattering perturbative calculations, arbitrary cou- 
pling constants gi, 52 > 93 for each type of vertices to keep track of perturbative calculations. As in the 
two-particle scattering case, we should take the limit gi ^ 1, i = 1,2,3 in the end. As we will show below, 
this is a crucial point, making the S-matrix factorization possible. Clearly, dealing with more vertices in this 
case is not an easy task already for the three-particle case even at the one-loop level. As the number of 



particles considered in the scattering increases, one has to expand the action ( 2.10 ) to higher order, which, as 
a result, makes the complexity of the perturbative calculations and the combinatorial analysis practically im- 
possible. Nevertheless, we make the first step in this direction and show in the next section the three-particle 
S-matrix factorization in the first order of perturbation. 

3.1 Diagrammatic calculations 

Assuming analyticity of the three-particle scattering amplitude in the coupling constants gi it can be written 
as: ^ 

(k|5|p) = (k|p) + (k|5|p) +{k\S\p) +••• (3.3) 

9 9^ 

where k = {ki, k2, k^}, and g is either of gi, i = 1,2, 3. The non-scattering term has the form: 

(p|k) = 3\{27TfSp[5{pi - ki)5{p2 - k2)6{p3 - ks)], (3.4) 
while the tree-level part yields: 

(k|5|p) =9ii2TTfSk,p{2giikik2+piP2)27r6ik3-p3)+g2ikik2+piP2)-2g3kipi}6E5P (3.5) 
9 

We have used here the symmetrization operator defined by [1]: 

'5p[/(p)] = |E/(^P) (3.6) 

A 

where the sum is taken over all possible permutations of (1,2,3), and the vector ^p = {pAi,PA2jPAs) is the 
corresponding set of momenta.^ We have also introduced 6E and 6P to denote the delta functions for the 
total energy and momentum conservation, i.e.. 



^In the case of the NLS model, the analyticity of the S-matrix can be seen from the N-particle coordinate Bethe ansatz 
calculations of Yang [64,65]. In the case of the LL model it is a natural assumption. 

^Without any loss of generahty, we choose the momenta p to be arranged in the order pi > P2 > Ps- 
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6E = s(j2{k^-p')^ 



(3.7) 
(3.8) 



The tree level Feynman graphs (connected and disconnected) in the first order are pictured in Fig. |2[ where 
the spatial derivatives are represented by the marks over the particle lines. We stress that each of these graphs 
denotes a sum over diagrams of the same topology but with different permutations of the external momenta. 



91 






+92 





Figure 2: The tree-level Feynman diagrams for the three-particle scattering in the Landau-Lifshitz model. 



The expression for the tree level (|3.5D defines the interaction vertex in the momentum representation, 
represented by Fig. 0. 




(k|5|p) 



Figure 3: Diagrammatic representation of the interaction vertex in momentum space. 

Using the interaction vertex in the momentum representation we are able to compute the second order 
scattering amplitude as depicted in Fig. ^ 




Figure 4: Diagrammatic representation of the second order scattering amplitude. 



The analytical expression corresponding to the Fig. H is given by: 
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(k|5|p) 



6 



n 

j=l L 



9z(27r)25fe,g{2sti(/ci/c2 + giQ2)2vr(5(fc3 - gs) + g2{kik2 + 91^2) 



+ g2{kik2 +P1P2) - 2<73A:iPi} 5 1^5^(9° 6 {^{qi-Pi)^ (3.9) 

The delta functions in ( |3.9| ) may be used to integrate over 93 and 93, and the integrals over q^ and 92 can 
be performed by choosing the contour closing in the lower half-plane. We also symmetrize over q, in order to 
obtain an analytical expression depending only on the external symmetrized momenta, 



(k|5|p) 



3i 



5E5PS. 



+ gi/2(k, q)/2(p, q) + ^glkipi ^k^p + 47r5ig2(/l(k, q)/2(p, q) + /2(k, q)/i (p, q))- 

- 87r5i53 (^Pi ^p/i(k,q) + /ci ^ /c/i(p, q)^ - 

- 25253 (pi^p/2(k,q)+A;i^A;/2(p,q))]} (3.10) 
where we introduced the following functions: 



dqidq2 



[I67r2g2/^(k,q)/i(p,q) + 



/i(x, q) = {xiX2 + qiq2)5{xi + X2 - qi - 92) + 



X1X2 + 91 



(E 



X- qi- q2 



Six3 - q2) + 



+ 



X1X2 + 92 - 91 - 92) (5(x3-9i) 



/2(x, q) = 3xiX2 + (91 + 92) ^ 



2; - 9l - 92 - 91^2 



(3.11) 
(3.12) 



To avoid cluttering we do not write the indices in the sum, and use the notation ^x = X^Xj. 

i=l 



It is easier to compute each term of (|3.10D separately. 



(a) Term proportional to gf: 

Most of the integrals are trivially calculated due to the delta functions in /i(x, q) , and the remaining 
delta functions corresponding to the total energy and momentum conservation allow us to write these 
terms as a sum of the following contributions: i) the finite part, which is a function of the external 
momenta, and ii) the term proportional to 6{ps — k^) . It is the latter term that may lead to divergencies, 
which, however, could be regularized in the same manner as in [27]. Indeed, we notice that this term 
can be written in the following form: 



-ikik2PiP2lo{ki, k2) + i{piP2 + kik2)h{ki, ^2) - ^-^2(^1,^:2), 
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Figure 5: Feynman diagram corresponding to the equation ( 3.13 ) 



where the one- loop integrals Iq, Ii and I2 were regularized and explicitely calculated in [27]. (We refer the 
reader to [27] for definitions and all calculational details). The finite result after this regularization has 
the form (where we keep here the integral representation to be dealt only in the end of the calculations) : 

gg^2 2 +P3{h + k2 -P3)]\P1P2 + h{Pl +P2- h)] 

^p0-p2_/c|-(pi+p2-/c3)2+ze 

2 2 f, [^-2piP2][^ + i(pi-p2)^-i(fcl-fc2)^-2fci/c2] ,„^ 

- 247r 51 / dq S{k3 - ps) (3.13) 

J - \[pi - P2)- s - le 

We have denoted s = ^ i^P^ — X^P^) and shifted the integration variable q ^ q + V2(^i + ^2)- 
The corresponding diagrams are presented in Fig. ^. Once again, we emphasize that each diagram 
corresponds to the sum over all diagrams with the same topology, but with distinct permutations 
over the external momenta and different placements of the derivatives. The first diagram in Fig. ^ 
corresponds to the first term in ( |3.13| ), and the second term is the finite part of the second diagram. 
The diagrams, corresponding to ^, but with explicit placement of the derivatives are depicted on Figs. 
^ and respectively.^ 

(b) Term proportional to g\: 

After some algebraic manipulations this term can be written in the form: 



2 f /2(k,q)/2(p,q) dqidq2 

^^7 T.P^ - Oi - <ll - {Y.P - 11 - 12^ + ie 

5i /■ , , [3piP2 - s] [3kik2 - s] 

dQidq2 2^2^ 7 — T — ' (3- 14) 

9i + 92 + 9192 -{qi + q2)2^p-s-ie 



The corresponding Feynman diagram is represented by Fig. 

(c) Term proportional to g^: 

Simple transformations lead to the following expression: 



. 2 f kipiY^kY^p dqidq2 



Ep° - Q1-Q2- (Ep- 91 - 92)^ + «e 



-2gl I k,p,EkEpdq,dq2 ^3^^^^ 

9i + 92 + 9192 -{qi + q2)2^p-s-te 



The corresponding Feynman diagram is represented by Fig. 



'We draw the remaining Feynman diagrams in the Appendix. 
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(d) Term proportional to gig2' 

One of the integrals in this term is trivially computed using the delta functions of /i(x, q), while 
the remaining integral after some algebraic manipulations and the shift of the integration variable 
qi ^ qi + y^ih + ^2) and 92 ^ 52 + V2(pi + P2), yields: 



4715152 / dqidq2 



/l(k,q)/2(p,q) + /2(k,q)/i(p,q) 

- qf - - (.Y.P- Qi - q^Y + «e 

6715152 / dq — ■ 

- + ^2)^ - ^3(^1 + k2)-s-ie 

« , , [p-iipi + P2) - P1P2 + s\[s - ?,kik2\ 
- 6715152 / dq — -( — - — V2 -( — - — ^ (3-16) 

2 - aVPi +P2) -P3[Pi +P2) - s-ie 



The corresponding Feynman diagram is represented by Fig. 10. 

(e) Term proportional to 5153: 

Once again the delta functions in /i(x, q) allow us to trivially compute one of the integrals and the 
resulting expression has the form: 



-8715153 / dqidq2 



Y.p° -qf-ql- iY.p-qi- 92)^ + ie 

kiY.^ [P3(Pi + P2) - P1P2 + s] 
q^ - \kP\ +P2Y -PsiPi +P2) - s-ie 



Ungigs / dq 



10 /J PiT,P[hih + k2)-kik2 + s] 
- 127r5i53 / dq- (3.17) 

q^ - j{ki + k2y - k3{ki + k2) - s - te 



The corresponding Feynman diagram is represented by Fig. |11 

(f) Term proportional to 5253: 

There are no delta functions present in this term and after some transformations we write it in the form: 



o , l'iEp/2(k,q) +A;iX;/c/2(p,q) 
25293 / ^F^-f^ 2 2 7^^ . dqidq2 



/ , , iPiEP + kiEk)s- 3kik2Pi E P - 3fcipip2 E k , . 

-5253 / aqiaq2 — -, ■ r^F^ : (3.i»j 

qi+qi + qiq2-[qi+q2)2^p- s-ie 



The corresponding Feynman diagram represented by Fig. 12 



Let us note, that in equations ( ^) , ( ^lej) , ( ^) and we have ignored the divergences 

and kept only the finite parts. Although the above formulas have been obtained for the general off-shell 
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case, putting the external momenta on shell, i.e., = pf and k!^ = kf, for i = 1,2,3 leads to considerable 
simplifications. This is due to the following identities: 



2{piP2 + PiP3 + P2P3), - (X^^) = -'2:{kik2 + kiks + k2k3) 



(3.19) 



+ P1P3 + P2P3 = kik2 + kik-i + k2k2, 



The remaining integrals are not difficult to compute. In fact, the integral that one has to compute in the 
second term of (|3.13|) is simply. 



dq 



2TTi 



q"^ - j{Pi - P2y - ie \pi-P2\ 



The remaining integrals in ( 3.16| ) and ( 3.17| ) are: 



(3.20) 



9^ ~ 3(^1 + ^2)^ - ^3(^1 + ^2) + kik2 + ^1^3 + ^2^:3 - 



dq 



9^ - jiPi +P2Y -PsiPi +P2) +P1P2 +P1P3+P2P3 - «e 



\ki - k2\ 

2m 
\Pi -P2I 



(3.21) 
(3.22) 



Finally the remaining double integral in ( p. 14 ), ( 3.15| ) e ( 3.18| ) is: 



/ 



dqidq2 



Q1+Q2 + - P{qi + 52) + iPiP2 + PiPs + P2P3) - if- 



7! 



^-ilogg2+l lini log(i?2-Q2) 



(3.23) 



where we denoted: 



= g [(Pl - P2? + (Pl - P3f + (P2 - P3?\ 



(3.24) 



Keeping only the finite part in the integral ( 3.23| ), and collecting the above expressions, we find that the 
second order scattering amplitude has the form: 



(ki^ip; 



9^ 2 
+ k\p\ 



[kik2+P3{ki + k2 -P3)][PiP2 + k^iPi +P2- k^)] kik2PiP2 



pl+pl-kl- {pi +P2- kzY + ie 
k2Pi \ in 



2m5{p3 - ks) 



+ kik2PiP2 



P1-P2 \ki-k2\ 
1 



+ 



^/3 



Pi -P2 

(ivr - logQ^) {{kip2 + k2Pi)g3{2gs - 3c/2) + g^kipi) 



+ 



+ 



Pi -P2 



+ 



^j-^ 51 (253 - 3ff2) + ^{iTT - logQ2)(253 - 352)' 



5E5P 
(3.25) 
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3.2 The S-matrix factorization 



Had we assumed the quantum integrability of the LL model, we could have obtained the S-matrix for three 
particle scattering by using the relation: 



{k\S\p) 



5(p)jf (p,k) 



(3.26) 



where 6^^\p,'k) = (p|k) = 3l{2iT)^Sp[6{pi — ki)6{p2 — k2)5{pz — ^3)]. This relation, analogous to ( p.lSD , 
underlines the fact that in the scattering process the particle annihilation and creation are not possible, and 
the set of momenta before and after the scattering is the same. This is not obvious from the calculations 



presented in the previous section. Indeed, although the non-scattering term clearly has the form ( 3.26 ), neither 
the tree level (^.5|), nor the second order ( 3.25 ) terms are manifestly of this form. Nevertheless, we will show 
that this is indeed the case, and one can reduce the expressions for the tree level (|3.5D and the second order 
( 3.25| ) contributions to the forms that are proportional to (5^^(p, k). 

Let us first focus on the tree level term ([3.5|). Expanding explicitly the symmetrization operator S^^p-, we 

find: 



(k|5|p) 



2igi 



P1V2 



+ 



+ 



<^?(p,k)+ 



LPI-P2 Pl-PS P2-P3_ 

+ 2i{2Tif {{pip2 + P1P3 + P2P3)(3(?2 - 2^3) - gsipl +pI+ pI)] 5E5P 



(3.27) 



We immediately notice that the two terms in above expression are quite different. Namely, the first term of 

(|3.27] ) is already in the form needed for the S-matrix to satisfy ( p.26| ), while the second term, on the other 

hand, does not satisfy this condition. Therefore, had we restricted ourselves only to this order, the condition 

(|3.26| ), necessary for integrability to hold, would have failed. Similarly, it is not difficult to see that the second 
^ f 3^ 

order scattering amplitude also contains terms that are not proportional to (P)l^)) such as the terms in 
the second and third square brackets in ( p. 251 ). 

Let us note, however, that if the quantum integrability, and as a consequence, the S-matrix factorization 
were true, the three-particle S-matrix would have the exact form: 



5(3) (p,k) 




+ o(r 



(3.28) 



where we have expanded in the last line the exact expression up to the second order in g. Comparing the 
equations (|3.27 ) and ( 3.2j ), we see that the unity in ( |3.28|) corresponds to the non-scattering part (p|k), 
while the first term in (|3.27l ) corresponds to the first order term of ( 3.28 ). We will now demonstrate, that 
the second term in ( |3.27| ) is such that it can be canceled out exactly with a term from the second order 
contribution ( 3.25[ ). 
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Figure 6: The Feynman diagrams for the first term in the equation (3.25) with all derivatives placed. 



A very similar expression for first term of ( 3.25| ) appears in the NLS model (see formula (13) of [1]). In 
the latter case, the use of the formula 

^ TiTT5{x)+F.V.(-] (3.29) 



xziziO \x ^ 

proved to be crucial to show the factorization (see also [53]). We, thus, expect that the usage of this formula 
should play the same role in our case. However, while the principal value part has no contribution for 
the NLS model, in the LL model it plays a non-trivial role and leads to the cancellations which make the 
factorization possible. We draw all the Feynman graphs for this term in Fig. ^, where only the sum over 
topologically similar diagrams with different permutations of the external momenta is implicit. Thus, there 
are four diagrams, which correspond to all possible placements of the derivatives. 

We remind that we have chosen the momenta p to be arranged in the order pi > P2 > Ps for the three- 
particle scattering to be possible. The principal value of the first term in ( |3.25| ) can be computed if we take 
ki / pj, i,j G {1, 2, 3} and use the identity 

1 i / 1 

+ 



kf + ^ -pf - {ki + kj - piY + ie 2 ki - kj \pi - kj - ie h - pi 



as well as take into account the delta functions corresponding to the energy and momentum conservation. 
Only the first and last diagrams on the R.H.S of the Fig. ^ will contribute to the principal value, while the 
other two vanish identically due to the momentum conservation. After some tedious transformations one 
finds: 



I = j^(~Pi -P2-P3+ P1P2 + pm + P2P3) 

(3.30) 

The delta function part of ( |3.29|) yields: 

S +pI- kl - {pi +P2- fca)) = W7-^ r [S{pi - k^) + 5{p2 - A^s)] (3.31) 

2{pi - P2) 

Substituting this result back into ( |3.25| ) we find: 



P.V. <^ cS, 



[kik2 + psiki +k2- Pz)\[PiP2 + ksipi + P2 - k^ 



kn 



(pi +P2- k-iY + ie 
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(k|S|p) 



%5E5P Sp,k \-96n^gl^i^^P^2ni [5{h - Pi) + 5(^3 - P2) + - ps) 
^ L Pi — P2 



+ hPi 

+ hk2PlP2 



ASir'^igigs ( ^^^^ + , , , ) + ^(^tt - logQ^) ((fcip2 + k2Pi)g3{2gz - 852) + glkiPi] 

P1-P2 \ki - A;2|y V3 



1 



+ 



1 



47r 



51(253 - 352) + ^(ivr - logQ2)(253 _ 352)' 

PI-P2 \kl-k2\J ^/3 



(3.32) 



To proceed, we explicitly expand the symmetrization operator Sk^p and obtain the following form for the 
second order scattering amplitude : 



+ 



(k|5|p) 



8ivr g^ (pi +P2+P3- P1P2 - P1P3 - P2P3)SESP+ 

2 



+ 2i^gl 



P1P2 



_^ PlP3 _^ P2P3 



iPl -P2 Pi- P3 

where we have separated the one loop contribution: 



P2 -P3 



6^^\p,k)-^n{g')6E6P 



(3.33) 



16 



-71" igi 



P1P2 ^ PlP3 ^ P2P3 



+ 



[(392 - 2g3){kik2 + kik3 + k2k3) - g3{ki + kl + kl)] + 

Pi -P2 Pi -P3 P2 -P3\ 

[(3^2 - 2g3){piP2 + PlP3 + P2P3) - g3{Pl +P2+ pI)] \ + 



kik2 ^ kiks ^ k2k3 



\ki-k2\ {ki-ksl 1^2-^3! 

+ — (^TT - logQ^) [(352 - '2g3){PlP2 + PlP3 + P2P3) -g3{Pl +P2+P3)] 



V3' 

[{3g2 - 2g3){kik2 + kiks + A;2A;3) - g3{kj + kl + h 



(3.34) 



We note that in the limit 91 = 52 = 93 = 1 the first term of ( 3.33| ) has the opposite sign of the second term 
in ( |3.27| ). Thus, the unwanted term indeed cancels out and the remaining terms combine into the form: 



(k|^|p) 



n=l 



P1P2 



+ 



PlP3 



+ 



P2P3 



Pi -P2 Pi- P3 P2 - P3 



3i. 

~2 



5f{p,\<i)-^n{g')5E5P + 0{g' 



(3.35) 



The first term in ( 3.35| ) is the three-particle scattering S-matrix in the second order, consistent with the 
factorization relation ( |3.28D . Thus, we conclude that to have the S-matrix factorization already in the lowest 
order, the terms from the higher order contribution should be taken into account to cancel out the unwanted 
terms. Diagrammatically, by making use of the relation ( ^.29 ) we have found that the principal value of the 
Feynman diagrams depicted in Fig. ^ cancels out the terms at first order that prevented the factorizability 
of the S-matrix at this order, while the contribution of the delta- function term of ( 3.2S| ), added to the term 
in Fig. 1^, gave the factorizable second order contribution to the S-matrix. Although we left out the TZ{g'^) 
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contribution in ( 3.35 ), this calculation provides us with a remarkable scheme in which the contributions of 
different orders cancel each other to yield S-matrix factorizability. 

It is not difficult to see that the diagrams that we had to sum up to obtain the three-particle S-matrix in 
the first order are of the order since neither the tree level, nor the diagrams in Fig. ^ contain any loops. 
Clearly, the cancellation at higher orders may happen for the diagrams of the same order in h. Even though 
we were able to prove the factorizability at the first order, the complexity of diagrammatic and combinatorial 
analysis make the computations of higher orders practically impossible to carry out, and another approach 
has to be employed. This subject is currently under investigation. 



4. Conclusion and discussion 

In this paper we have shown the three-particle S-matrix factorization at the first order for the LL model. The 
calculations we have presented show the non-trivial mechanism behind the factorization even in the lowest 
order of perturbation. In the process we have also shown the absence of particle annihilation and creation 
in the scattering processes, and that the set of momenta before and after the scattering, correspondingly, is 
unchanged. Our consideration was motivated by the difficulties to analyze the quantum integrability of the 
LL model in the framework of the quantum inverse scattering method. The highly singular nature of the LL 
model does not allow the use of the trace identities and, as a result, the set of local operators in involution 
is in general hard to derive. In fact, only the action of the quantum Hamiltonian on the two-particle sector 
is known. In addition, the standard two-particle state underlining the clear meaning of the scattering 

process, is not valid for the LL model quantized in terms of the fields (p{x) ( |2.9D , namely, the particles 
created by the operator ip~^{x) are not the Bethe particles. As the S-matrix factorization, which expresses 
quantum integrability, relies on the presence of an infinite set of commuting local operators, it is important 
to independently verify quantum integrability using field theoretic methods. Let us note, that unlike the LL 
model, these subtleties are absent in the classically equivalent NLS model, and the quantum inverse scattering 
method is easily shown to be consistent with the field theoretic calculations. 

The next obvious problem is generalizing the above construction to all orders of perturbation, as well 
as showing the factorization for the A^-particle scattering amplitude. This is not an easy task, as even in 
the lowest order the difficulties of diagrammatic and combinatorial nature are quite complex. The essential 
difference with the NLS model is the new vertices that appear in each order for the LL model. As we have 
shown, the factorization mechanism is such that the higher order terms cancel the unwanted terms via the use 
of the formula ( 3.29| ). Clearly, at higher orders the analysis will be too complicated to be carried out. Thus, an 



alternative framework is needed to perform the diagrammatic calculations. The latter have to be also carried 
out to understand the effects of renormalization. The reason is twofold. First, as we have shown, the coupling 
constants gi,i = 1, 2, 3, which we had introduced by hand to keep track of the perturbative order, had to be set 
to unity for the three-particle S-matrix to be factorizable. Although we have done so without properly taking 
care of the renormalization, it is important to carry out all the necessary calculations to show the validity of 
this assumption. Secondly, it is clear that the singularities arising in the quantum inverse scattering method, 
when deriving the local commuting operators, should exhibit themselves in the renormalization of the coupling 
constants and fields. Thus, it would be interesting to establish a direct map between the two methods. We 
stress the importance of this program in the context of the AdS/CFT correspondence, since the difficulties 
present in the quantization of LL model will also be present in the higher sectors of the AdS^ x string. 
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Appendix: Feynman diagrams 

In this appendix we draw the remaining Feynman graphs for the three-particle second order scattering ampli- 
tude with the explicit placement of derivatives. We note, however, that the sum over the diagrams of the same 
topology, but with different permutations of the external momenta, is assumed. The two diagrams depicted 
in Fig. m, which correspond, respectively, to each of the finite terms of the equation ( |3.13| ), are represented, 
with all derivatives placed, by Figs. |6| and ^. 

The following list of figures contains the complete Feynman diagrams corresponding to the equations 

dull)- 



Figure 7: The Feynman diagrams with the correct placement of derivatives corresponding to the second Feynman 
graph pictured on Fig. ||. 
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+ 



(3.15) 



+ Time Reversed Diagrams 

Figure 10: The Feynman diagrams proportional to (71(72, correspon(iing to equation ( 3.16|) . 
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